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EFFECTIVE THERMAL CONDUCTIVITY DETERMINATION 

FOR LOW-DENSITY INSULATING MATERIALS 

By S. D. Williams* and Donald M. Curry 
Lyndon B. Johnson Space Center 

SUMMARY 


A nonlinear least squares technique has been used to determine effective 
thermal conductivity values using experimental data for a low-density insu- 
lating material. A comparison of the analytically predicted conductivity 
with values obtained from standard guarded-hot-plate tests indicates that the 
predicted values can be used with confidence in performing thermal analyses. 

An assessment of the predicted thermal conductivity, from a statistical view- 
point, and of its effect on the thermal response of the material was performed 
and is discussed herein. Results of this investigation demonstrate that ther- 
mal conductivity values can be routinely predicted using analytical methods 
rather than being obtained through expensive, time-consuming standard guarded- 
hot-plate tests. 


INTRODUCTION 


The design and development of a reusable thermal protection system (TPS) 
for the Space Shuttle Orbiter required a detailed investigation of various 
material classes and configurations, and a knowledge of the thermophysical 
properties of these materials. Since direct measurement of effective thermal 
conductivity is both expensive and subject to inaccuracy, guarded-hot-plate 
(6HP) data may be supplemented with effective conductivity inferred from in- 
depth thermocouple (TC) temperature histories obtained from thermal evalua- 
tion tests (plasma arc jet/radiant). Using analytical methods, properties 
can be routinely derived from these types of tests rather than being obtained 
from expensive, time-consuming standard tests. 

Of special interest are analytical methods that permit the evaluation of 
thermal response data obtained from a single experiment. Numerical tech- 
niques have been developed to determine thermophysical properties and surface 
temperature response from in-depth temperature measurements. In particular, 
Curry and Williams (ref. 1) have shown that meaningful effective thermal 
conductivity can be obtained by use of nonlinear least squares techniques 
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from thermal tests not specifically designed to measure the thermal conduc- 
tivity. Banas and Cunnington (ref. 2) have discussed laboratory measurements 
of effective thermal conductivity, fiber diameter, and extinction coefficient 
for the all-silica, low-density, rigid insulation material selected for the 
Orb iter TPS. 

The application of the nonlinear least squares method in calculating 
effective thermal conductivity values using transient thermal response data 
is discussed in this paper. A comparison of the analytically predicted 
conductivity with values obtained from a standardized GHP apparatus (ref. 3), 
is presented. An assessment of the predicted effective thermal conductivity, 
from a statistical viewpoint, and of its effect on the thermal response of 
the material is also discussed. In addition, the results obtained from using 
a Space Shuttle mission profile response are given. These simulations dem- 
onstrate that reliable results can be obtained using derived conductivity 
values rather than GHP data. 

As an aid to the reader, where necessary the original units of measure 
have been converted to the equivalent value in the Systeme International 
d'Unites (SI). The SI units are written first, and the original units are 
written parenthetically thereafter. 


SYMBOLS 


ao,ai,a 2 ,a 3 coefficients in conductivity polynomial (eqs. (2) and (5)) 

dYi error between predicted and measured thermocouple data 

(eq. (8)) 

kc conductivity due to conduction and convection (eq. (1)) 

kg effective thermal conductivity (eq. (1)) 

kp conductivity due to radiation (eq. (1)) 

Z degree of polynomial fit (eq. (8)) 

m number of time measurements (eq. (3)) 

n number of thermocouples (eq. (3)) 

N number of baseline conductivity values (eq. (9)) 

p number of measurements taken with a thermocouple (eq. (8)) 

Q-j^ input heat (fig. 2) 

r level of confidence (eq. (7)) 

S sum of squares (eq. (3)) 
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T temperature (eq. (2)) 

T(t) temperature function (fig. 2) 

Tij predicted temperature at i-th thermocouple and j-th time 

measurement (eq. (3)) 

T-jj* measured temperature at i-th thermocouple and j-th time 

measurement (eq. (3)) 

X depth of test article (fig. 2) 

Akq difference between baseline conductivity and predicted 

effective conductivity (eq. (9)) 

Ak statistical mean of Akq (eq. (9)) 

standard error of estimate of conductivity (eq. (5)) 

standard deviation of conductivity between baseline and 
predicted effective conductivity (eq. (9)) 

standard error of estimate of temperature (eq. (4)) 

standard error of estimate of temperature for thermocouple 
data (eq. (8)) 

Subscripts 

i thermocouple measurement 

j time measurement 

THEORETICAL FORMULATION 


To evaluate an effective method for determining thermal property values, 
it is necessary to establish the fitting equations that are to be used. The 
purposes of the curve-fitting equations are to summarize a large quantity of 
data, to establish interpolation equations, and to confirm the analytical 
model . 

Since it may be difficult to obtain explicit comparisons with the fit- 
ting equations, criteria should be established to satisfy the analyst that 
the results are reliable. The following criteria were used to ensure a good 
fitting method. 

1. Use all of the relevant data in estimating each constant. 

2. Have reasonable economy in the number of constants required. 
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3. Provide some estimate of the uncontrolled error in the prediction. 

4. Make it possible to find regions of systematic deviations or func- 
tional bias from the equation, if any exist. 

5. Show whether the conclusions are unduly sensitive to the results of 
a small number of runs, perhaps even to one run. 

6. Give some idea of how well the final equation can be expected to 
predict the response, both in the overall sense and for important sets of 
conditions inside the region covered by the data. 

The material investigated in this study is a low-density (144.2 kg/m3), 
rigid, fibrous insulation with a high porosity (94 percent void volume) and 
with transparency that permits heat transfer to occur by conduction, con- 
vection, and radiation. The effective thermal conductivity kg, discussed in 
this paper can be expressed as 


k 


e 





( 1 ) 


where k^ is the conductivity due to conduction and convection, and kR is 
the apparent conductivity due to internal radiation transfer. In this inves- 
tigation, it is assumed that the effective thermal conductivity can be 
expressed as a polynomial in temperature 


'e = 


ag + a^T + ^ 2 '^ 


asT' 


( 2 ) 


where the values of the unknown coefficient ai are to be determined for a 
minimal difference between the measured temperatures and the predicted 
temperatures. If Tij represents the predicted temperature at the i-th 
thermocouple and the j-th time measurement, the least squares problem is to 
minimize 


m 

S ^ag,a^ »a2’®3) ^ S 

j=l 


n 



T. . - T. .* 



(3) 


where S is the sum of squares, m is the number of time measurements, n 
is the number of thermocouples, and T-jj* is the measured temperature. Ref- 
erence 1 provides the justification for selecting the model in this form and 
demonstrates that the most efficient technique in solving for the unknown 
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coefficients is one developed by Peckham (ref. 4). It is obvious that the 
first criterion is satisfied by the selection of thermocouple locations and 
the choice of times during the test at which data will be selected for fit- - 
ting. Since only four coefficients are used, the second criterion is also 
satisfied. 

The third criterion can be satisfied by examining the converged least 
squares fit. A statistical measurement can be made to determine the accuracy 
of the converged effective conductivity polynomial. This analysis is given 
in the form of the standard error of estimates of temperature o-r and con- 
ductivity from the least squares fit of the data. The standard error 
of temperature is given by 
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where the number 4 accounts for the number of coefficients used to fit the 
data. To determine the standard error of conductivity from and the rep- 
resentation for conductivity, one has 



9k 




e 

9T 

II 

Q 

r+ 

a«| + 2d2T + 


It has been found that the values of alone offer little meaning to the 
analyst and cannot be effectively used to satisfy the fourth criterion. How- 
ever, if one uses °k as a measure of the accuracy of kg, the information 
is useful. A convenient method is to use °k ^ percent error of the ef- 
fective conductivity; thus. 


Percent error = ^ x TOO (6) 


This method can also be used to establish confidence limits on the accuracy 
of the data fit. Since there is no guarantee of the least squares fit being 
normally distributed, one must use the more conservative Chebyshev inequality 
(ref. 5) 
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(7) 


±ra|^ ^ ^ 100 


where r is the level of confidence. In terms of the percent error 
expressed in equation (6), this expression means that a 3o confidence level 
(90 percent) on the accuracy can be related in terms of conductivity. Thus, 
if one has a 3-percent error in conductivity, he can be reasonably assured 
that a 9-percent error will bound 90 percent of the cases investigated. 


Test Data Reduction 

In the least squares procedure, it is assumed that the independent vari- 
ables are measured without error. However, the data resulting from thermal 
simulation tests always contain a certain amount of noise and usually have a 
few "wild" points (i.e., data values that are so far removed from the trend 
of the other data that they are obviously invalid). To the analyst, these 
irregularities present no problem since the data are simply plotted at inter- 
vals much greater than the sample rate, and obvious wild points are discarded. 
Synthesis of the thermal response data requires that the computer program 
accept raw thermocouple data as input, usually from a magnetic tape. Diffi- 
culties arise, however, when the raw thermocouple data are used since, during 
the course of the repetitive data reduction and thermal calculations, even 
small errors can be magnified to make the output meaningless. 

An algorithm has been developed to smooth the thermocouple data by 
fitting various least squares polynomials on overlapping data subsets. An 
attempt to fit all of the temperature data with a single polynomial would 
lead to large errors since a single polynomial cannot be found that repre- 
sents all of the data. In addition, the fit at the end points would be ex- 
tremely poor. 

Exhaustive numerical testing has shown that the best results can be 
obtained by using a seventh-degree polynomial on each segment for the least 
squares fit. When all data are collected, the standard error of temperature 
Oy can be computed from 
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i=l 


where dY-j is the error computed, p is the number of thermocouple data 
measurements, and I is the degree of the polynomial fit. The standard 
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error of temperature Oy can be used to accept or reject thermocouple data 
and thereby to ensure use of only those data that will produce meaningful 
results. 


Thermocouple Location Errors 

It is often difficult to physically locate accurately the position of 
the thermocouples; however, in predicting thermal conductivity by means of 
the mathematical model, these locations are assumed to be exact. After con- 
verging to the least squares approximation, the nonlinear least squares 
(NONLIN) program attempts to relocate the thermocouples in order to reduce 
the least squares errors. Frequently, this procedure alters the conductivity 
values such that they become more representative of the actual conductivity. 
In addition, by analyzing the results from the relocated thermocouples, the 
fifth criterion can be satisfied in an economical manner. 


Prediction Technique Validation 

The validation of the predicted effective thermal conductivity values 
can be established by several procedures. The most obvious method is to use 
the results to predict in-depth thermal responses and compare them with meas- 
ured data. The next procedure is to use the conductivity values obtained 
from different constant-pressure tests to predict the in-depth response of 
the material subjected to variable heating and pressure conditions. The 
final validation consists of using the conductivity values obtained from the 
test specimen to predict the in-depth response of additional test specimens. 
These three validating procedures will satisfy the sixth criterion. 


EXPERIMENTAL/ANALYTICAL VERIFICATION 


A series of thermal tests using the Space Shuttle Orbiter reusable 
surface insulation (RSI) material was performed to validate the design ther- 
mal properties. 


Test Conditions 

The tests consisted of exposing the RSI TPS test article to steady-state 
temperature conditions of 1367 K at various constant pressures over a range 
from 0.10 to 101 kN/m^ (0.76 to 760 torr). In addition, the test articles 
were exposed to transient trajectory temperature and flight pressure profile 
conditions. The test chamber was evacuated to the desired test pressure 
condition for a minimum of 1 hour before initiation of the heating cycle. 

Each test cycle was initiated when all thermocouples were essentially at the 
same temperature. Data were taken for a minimum of 1 hour after initiation 
of the heating cycle. Heater operating conditions, environmental test 
conditions, and model temperature response data were recorded on magnetic 
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tape on an analog-to-digital recording system. The data on the tapes were 
reduced to engineering units by use of a Univac 1100 series computer system. 


Test Article Configuration 

A cross-sectional view of the test article/test fixture is shown in 
figure 1. Two RSI tiles (15.24 by 15.24 by 11.684 cm), which were fabricated 
using baseline silica fibers and were coated on the top surfaces and side- 
walls with a reaction-cured glass (RCG) formulation, were provided by the 
Lockheed Missiles and Space Company. Each tile contained a 3.81-cm-diameter 
plug that was instrumented with six 127-micrometer (5 mil) platinum/platinum/ 
13-percent-rhodium thermocouples in the center of the tile. 


EMBEDDED 




Figure 2.- Thermal model used 
to simulate test article 
configuration. 


Effective Thermal Conductivity 
Predictions 

The NONLIN program was used to cal- 
culate effective thermal conductivity 
values using the thermal response data 
acquired from the constant-pressure/ 
radiant-heating tests. The results of 
each radiant test were analyzed using 
a one-dimensional thermal mathematical 
model (fig. 2) developed to simulate 
the test article configuration. Since 
it is difficiult to locate the exact 
thermocouple position, the NONLIN program, 
after converging to the least squares 
approximation, assesses positional error 
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effects and then repositions the thermocouples analytically to reduce 
the cTj value. 

A composite plot of the predicted effective thermal conductivity as a 
function of temperature (with pressures of 101, 10, 1.0, and 0.10 kN/m^ (760, 
76, 7.6, and 0.76 torr) ) using the steady-state test data with nominal and 
relocated-thermocouple locations for specimen RB37 can be seen in figures 
3(a) and 3(b), respectively. A comparison with baseline thermal conductivity 
values is also indicated in these figures. These baseline values were ob- 
tained using a GHP apparatus with earlier material samples. The effect of 
using the relocated-thermocouple locations, in this case, was to reduce the 
values of by 50 percent (table I) and to raise the conductivity values 
slightly. A comparison of the calculated thermocouple depths with the nomi- 
nal values can be seen in table II. By relocating the thermocouples, the 
percent error based on also was reduced by more than 50 percent from a 
typical value of 5 percent to 2 percent. (See figs. 4(a) and 4(b).) This 
result indicates that a 3a confidence level of less than 6 percent error 
can be anticipated by using relocated data values. It is interesting to ob- 
serve that the percent error for lOl-kN/m^ (760 torr) and l.O-kN/m^ (7.6 torr) 
data is sufficiently large in the same region (less than 450 K) that an analyst 
would interpret the conductivity data as being subject to error. In addition, 
the values for are sufficiently large to make the data suspect. In gen- 
eral, the predicted conductivity values are approximately 18 percent higher 
than the GHP data^. However, the predicted conductivity values for pressures 
of 101 and 1.0 kN/m^ (760 and 7.6 torr) and temperatures less than 450 K are 
less than the GHP values. This reversal in calculated conductivity is not 
consistent with the other results but is explained by the larger deviation 
between measured and predicted temperature for thermocouples 5 and 6 at 101 
and 1.0 kN/m2 (760 and 7.6 torr). The results for thermocouples 5 and 6 are 
also consistent with the error deviations shown in table III between the raw 
test data and the least^squares polynomials. Table IV shows the standard devi- 
ation of conductivity between the baseline (GHP) and predicted effective 
thermal conductivity using both nominal and relocated-thermocouple depths. 


The standard deviation is 

defined as 
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Figure 3.- Comparison of predicted 
thermal conductivity for specimen 
RB37 (curves) with baseline con- 
ductivity values (data points), 
(overleaf) 


iThe GHP tolerances on the RSI 
thermal conductivity are ±18 percent. 
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Baseline (GHP) 
O101kN/m2 (760torr) 
□ 10 kN/m2 ( 76 torr) 
O 1.0 kN/iTi2 (7.6 torr) 
A .10 kN/m2 (.76 torr) 


Predicted (NONLIN) 

101 kN/m2 (760 torr), cry = 16.32 K 

10kN/m2( 76 torr), try = 17.33 K 

1.0 kN/m2 (7.6 torr), aj = 19.67 K 

.10kN/fTi2 (.76torr), <ry = 22.36 K 



Temperature, K 


(a) Nominal thermocouple locations. 


Baseline (GHP) 

O 101 kN/m2 (760 torr) 
□ 10 kN/m2 ( 76 torr) 

O 1.0 kN/m2 (7.6 torr) 
A .10 kN/m2 (.76 torr) 


Predicted (NONLIN) 

101 kN/m2 (760 torr), o-y = 8.06 K 

10kN/m2( 76 torr), try = 8.12 K 

1.0 kN/m2 (7.6 torr), <ry = 5.88 K 

.10 kN/m^ ( .76 torr), <ry = 9.52 K 



(b) Relocated- thermocouple locations 



TABLE I.- A COMPARISON OF THE STANDARD ERROR OF ESTIMATE 
OF TEMPERATURE FOR NOMINAL AND RELOCATED-THERMOCOUPLE LOCATIONS 


Pressure, kN/m2 Oj, K 

(torr) 




Nominal TC 

Relocated 

101.32 

(760.00) 

16.32 

8.06 

10.13 

(76.00) 

17.33 

8.12 

1.01 

(7.60) 

19.67 

5.88 

.10 

(.76) 

22.36 

9.52 


TABLE II.- A COMPARISON OF NOMINAL AND RELOCATED-THERMOCOUPLE DEPTHS 


Specimen Pressure, kN/m2 Type Thermocouple depths, cm 

(torr) (a) 






TCi 

TC 2 

TC3 

TC 4 

TC 5 

TCe 

RB37 

101.32 

(760.00) 

R 

0.0381 

0.376403 

0.997356 

1.636344 

2.747848 

5.456098 

RB37 

10.13 

(76.00) 

R 

.0381 

.378917 

1.002817 

1.630782 

2.699004 

5.473497 

RB37 

1.01 

(7.60) 

R 

.0381 

.389966 

1.013130 

1.611935 

2.710358 

5.510886 

RB37 

.10 

(.76) 

R 

.0381 

.391744 

1.016813 

1.634566 

2.668194 

5.223713 

RB37 

- 

- 

N 

.0381 

b.4572 

C.9398 

bl.6256 

^2.7686 

d5.3086 

B5892 

- 

- 

N 

.0381 

b.4572 

dl.5494 

^3.8608 

^6.5786 

C9.0424 


^R = relocated-thermocouple values; N = nominal values. 


*^Plus or minus 0.03 centimeter. 
*^Plus or minus 0.09 centimeter. 
*^Plus or minus 0.05 centimeter. 
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Conductivity percent error 


100 



200 400 600 800 1000 1200 

Temperature, K 


(a) Nominal thermocouple locations. 



200 400 600 800 1000 1200 

Temperature, K 


(b) Relocated- thermocouple locations. 
Figure 4.- Conductivity percent error based on a. 


TABLE III.- A STANDARD ERROR OF TEMPERATURE FOR LEAST SQUARES FITTING 
OF THERMOCOUPLE DATA FOR SPECIMEN RB37 


Pressure, kN/m2 Standard error of temperature, K 


(torr) 



I J 

TCi 

TC 2 

TCs 

TC 4 

TC 5 

TC6 

101.32 

(760.00) 

0.646 

3.335 

1.950 

2.422 

11.213 

5.878 

10.13 

(76.00) 

.343 

2.749 

1.764 

2.034 

9.966 

5.544 

1.01 

(7.60) 

.542 

3.763 

3.151 

3.151 

13.202 

7.197 

.10 

(.76) 

.511 

3.231 

2.058 

2.448 

12.723 

6.182 


TABLE IV.- 

STANDARD 

DEVIATION 

OF THERMAL 

CONDUCTIVITY 



Pressure, kN/m2 
(torr) 


Oj^, W/m*K 


Nominal TC 


Relocated TC 


101.32 

(760.00) 

0.0315 

0.0330 

10.13 

(76.00) 

.0182 

.0206 

1.01 

(7.60) 

.0192 

.0190 

.10 

(.76) 

.0150 

.0169 


where N is the number of baseline conductivity values, Akq is the dif- 
ference between 6HP and predicted conductivity values, and is the sta- 
tistical mean of Akq. Although the values obtained by using the 
nominal data are slightly less than those obtained by means of the relocated 
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data, the temperature errors are greater. The differences between the cf. 
values for the nominal and relocated data are not significant from a sta-^ 
tistical viewpoint; however, the differences in oj values are significant. 

Figures 5(a) to 5(d) present a comparison of measured and predicted in- 
depth temperatures for constant pressures of 0.10, 1.0, 10, and 101 kN/m2 
(0.76, 7.6, 76 and 760 torr), respectively, using nominal thermocouple 
locations and the corresponding predicted conductivity values for specimen 
RB37. Figures 6(a) to 6(d) present the same comparison of temperatures using 
relocated-thermocouple locations and the corresponding predicted conductivity 
values for specimen RB37. The predictions using the relocated data are in 
closer agreement with the measured data, which confirms the analysis based on 
both Oj and 0 |^. 

Because of the errors indicated in the predicted conductivity values 
(figs. 4(a) and 4(b)), engineering judgment was used to "fair" the con- 
ductivity values at pressures of 101 and 1.0 kN/m2 (760 and 7.6 torr). The 
same thermal response predictions were made as previously discussed, and the 
results agreed well with those presented in figures 5 and 6. This result 
confirmed that the data in this region were not significant in the thermal 
conductivity estimation. 

Thermal response predictions for the mission profile (variable surface 
temperature and pressure) using the predicted effective thermal conductivity 
and the baseline conductivity (GHP) for specimen RB37 are shown in figures 7(a) 
and 7(b), respectively. As can be seen, a better comparison between measured 
and predicted in-depth temperatures is obtained using the analytically deter- 
mined conductivity. The predicted thermal response for the last two thermo- 
couples (5 and 6) was lower than anticipated on the basis of the results 
of the constant-temperature/constant-pressure tests. This result is attrib- 
uted to the noise present for the first 600 seconds on all tests for these 
two thermocouples. 

Thermal response predictions for the mission profile using the predicted 
effective thermal conductivity and the baseline conductivity (GHP) for the 
second test article (specimen B5892) are shown in figures 8(a) and 8(b), re- 
spectively. The nominal thermocouple depths for specimen B5892 are given in 
table II. As before, a better comparison between measured and predicted in- 
depth temperatures is obtained using the analytically determined conductivity. 

These mission simulation tests demonstrate the uniqueness of the inferred 
effective conductivity values obtained from constant-surface-temperature/ 
constant-pressure tests. In addition, these comparisons demonstrate that 
more reliable results for low-density, high-porosity, and semitransparent 
insulative materials can be obtained by using analytically calculated effective 
conductivity values than can be obtained by using GHP data. The reasons 
for these better predictions using effective conductivity appears to be due 
to several factors. 

1. The GHP values are obtained from steady-state conditions and use a 
mean temperature in calculations for conductivity, whereas the analytical 
methods depend on local gradients in a transient environment. 
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(a) 0.10 kH/\n (0.76 torr). 


Measured 
O TC 1 
□ TC 2 
O TC 3 
A TC 4 
V TC 5 
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X = 0.0381 cm 
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(b) 1.0 kN/m^ (7.6 torr). 


Figure 5.- Comparison of temperature response using predicted thermal conduc- 
tivity from specimen RB37 tested at various pressures for nominal thermo- 
couple (TC) locations at depth X (curves) with measured thermocouple 
values (data points). 
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Temperature, K Temperature, 



Time, sec 

(c) 10 kN/m^ (76 torr). 


Measured 

Predicted (NONLIN) 


O TC 1 

X = 0.0381 

cm 

□ TC 2 

X = 0.4572 

cm 

O TC 3 

X = 0.9398 

cm 

A TC 4 

X = 1.6256 

cm 

V TC 5 

X= 2.7686 

cm 

TC 6 

X= 5.3086 

cm 





Measured 

Predicted (NONLIN) 


0 TC 1 

X = 0.038100 

cm 

□ TC2 

X = 0.391744 

cm 

0 TC 3 

X = 1.016813 

cm 

A TC 4 

X = 1.634566 

cm 

V TC 5 

X = 2.668194 

cm 

TC6 

X= 5.225713 

cm 



(a) 0.10 kN/m2 (0.76 torr). 


Measured Predicted (NONLIN) 

O TC 1 0.038100 cm 

□ TC 2 X - 0.389966 cm 

2 TC3 X^ 1.013130 cm 

TC 4 X ^ 1.611935 cm 

V TC 5 X - 2.710358 cm 

Tfr TC 6 X ^ 5.510886 cm 



(b) 1.0 kN/m2 ( 7.6 torr). 


Figure 6.- Comparison of temperature response using predicted thermal conduc 
tivity from specimen RB37 tested at various pressures for relocated- 
thermocouple locations at depth X (curves) with measured thermocouple 
values (data points). 
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Measured 

Predicted (NONUN) 

O TC 1 

X= 0.038100 cm 

□ TC 2 

X = 0.378917 cm 

O TC 3 

X= 1.002817 cm 

A TC 4 

X= 1.630782 cm 

y TC 5 

X = 2.699004 cm 

'k TC 6 

X = 5.473497cm 



Time, sec 


(c) 10 kN/m2 (76 torr). 


Measured Predicted (NONLIN) 

O TC 1 X = 0.038100 cm 

a TC 2 X = 0.376403 cm 

^TC 3 X = 0.997356 cm 

^TC 4 X = 1.636344 cm 

^ TC 5 X = 2.747848 cm 

TC 6 — — — X = 5.456098 cm 



(d) 101 kN/m2 (760 torr). 

Figure 6.- Concluded. 
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(a) Relocated-thermocouple thermal conductivity. 


Measured 

Predicted (NONLIN) 

O TC 1 

X - 0.0381 cm 

□ TC 2 

X = 0.4572 cm 

O TC 3 

X = 0.9398 cm 

A TC 4 

X = 1.62 56 cm 

V TC 5 

X = 2.7686 cm 

It TC 6 

X= 5.3086 cm 



Figure 7.- Comparison of mission profile temperature response at nominal 
thermocouple locations (curves) with measured values (data points) for 
specimen RB37. 


19 



1800 


1600 

1400 — 

1200 - 
oT 

f 1000 - 

Q. 

K 800 
600 
400 
200 


\ 


'•,0 . 




O' 


^ — O’ — "O • ~ O" 

,0^ r ‘ 

n 
/ / 

/ C3 / 

' / 

o' ^ 

--<1 1 ^— — * « r — 

-_1 i J . ^ u 

600 800 1000 1200 1400 1600 


O" * 
- " 


' ■■ ^ 


A' 


■ t: 

- 




- '(r 

. I 


0 200 400 

Time, sec 

(a) Relocated- thermocouple thermal conductivity. 


Measured 

Predicted (NON LIN) 

O TC 1 

X = 0.0381 cm 

□ TC 2 

X = 0.4572 cm 

O TC 3 

x = 1.5494 cm 

A TC 4 

— X = 3.8608 cm 

y TC 5 

X = 6.5786 cm 

« TC 6 

X = 9.0424 cm 



Figure 8.- Comparison of mission profile temperature response at nominal 
thermocouple locations (curves) with measured values (data points) for 
specimen B5892. 
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2. In-depth radiation effects are not fully accounted for under GHP 
test conditions but are closely approximated in the analytical model using 
transient data. 

3. The thermal conductivity values determined either through GHP or by 
analytical calculations are not exact, but the uncertainty in GHP values is 
probably greater than or equal to the uncertainties associated with the 
analytical calculations. 

4. The analytical calculations represent the best predictions over the 
entire temperature range at all times as opposed to isolated discrete 
measurements . 


CONCLUDING REMARKS 


That nonlinear least squares can be used with confidence to determine 
effective thermal conductivity from dynamic thermal performance data has been 
demonstrated. In addition, a method was presented for assessing the relative 
error associated with the predicted conductivity values. 

The thermal conductivity for a low-density, high-porosity insulation 
material was predicted from four constant-temperature/constant-pressure tests. 
It was demonstrated that the standard error of temperature and the percent 
error of thermal conductivity could be reduced by 50 percent through analyt- 
ical relocation of the thermocouple positions in the thermal model. Whereas 
the predicted data were 18 percent higher than guarded-hot-plate data, the 
error analysis indicated that there was less than a 6-percent relative error 
for a 3o confidence level in the analytically predicted data. 

An analysis of the thermal response for a simulated Space Shuttle mis- 
sion profile was made using both predicted and guarded-hot-plate thermal 
conductivity values for two separate test specimens. These tests demon- 
strated that more reliable results can be obtained using analytically deter- 
mined thermal conductivities. 

This report also points out the differences between dynamic (transient) 
and static (guarded-hot-plate) determination of effective thermal conductiv- 
ity where the means of heat transfer in the insulation material is by a 
combination of conduction, convection, and radiation. 

Lyndon B. Johnson Space Center 

National Aeronautics and Space Administration 
Houston, Texas, October 7, 1977 
986-15-31-04-72 


21 



REFERENCES 


1. Curry, D. M.; and wmiams, S. D: Nonlinear Least Squares > An Aid to 

Thermal Property Determination. AIAA J., vol . 11, no. 5, May 1973, 
pp. 670-674. 

2. Banas, R. P.; and Cunnington, G. R., Jr.: Determination of Effective 

Thermal Conductivity for the Space Shuttle Orbiter's Reusable Surface 
Insulation (RSI). AIAA Paper 74-730, July 1974. 

3. Standard Test Method for Steady-State Thermal Transmission Properties by 

Means of the Guarded Hot Plate. American Society for Testing and 
Materials Standard ASTM C 177-71, 1971. 

4. Peckham, G.: A New Method for Minimising a Sum of Squares Without 

Calculating Gradients. Computer J., vol. 13, no. 4, Nov. 1970, 
pp. 418-420. 

5. Winkler, Robert L.; and Hays, William L.: Statistics: Probability, 

Inference, and Decision. Second ed.. Holt, Rinehart and Winston 
(New York), 1975, pp. 169-172. 


22 



1. Report No. 2. Government Accession No. 

NASA TP- 1155 

4. Title and Subtitle 

EFFECTIVE THERMAL CONDUCTIVITY DETERMINATION 
FOR LOW-DENSITY INSULATING MATERIALS 

7. Author(s) 

S. D. Williams, Lockheed Electronics Company, and Donald 
M. Curry, Johnson Space Center 

9. Performing Organization Name and Address 

Lyndon B. Johnson Space Center 
Houston, Texas 77058 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

15. Supplementary Notes 


3. Recipient's Catalog No. 


5. Report Date 

February 1978 

6. Performing Organization Code 

JSC-13119 

8. Performing Organization Report No. 

S-481 

10. Work Unit No. 

986-15-31-04-72 

11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Paper 

14. Sponsoring Agency Code 


16. Abstract 


That nonlinear least squares can be used to determine effective thermal conductivity is 
demonstrated, and a method for assessing the relative error associated with these pre- 
dicted values is provided. The differences between dynamic and static (guarded hot plate) 
determination of effective thermal conductivity of low-density materials that transfer 
heat by a combination of conduction, convection, and radiation are discussed. 


17. Key Words {Suggested by Author(s) ) 

Heat transfer 
Thermocouples 
Standard deviation 
Mathematical models 
Statistical analysis 
Nonlinear least square s 

19. Security Oassif. (of this report) 

Unclassified 


18. Distribution Statement 


STAR Subject Category: 

34 (Fluid Mechanics and Heat Transfer) 


20. Security Classif. (of this page) 

Unclassified 


21. No. of Pages 

24 


22. Price* 

$3.50 


For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA-Langley, 1977 



National Aeronautics and 
Space Administration 

Washington, D.C. 

20546 

Official Business 

Penalty for Private Use, $300 


THIRD-CLASS BULK RATE 


Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 



\ 


IW\S«\ 


6 1 10, D, 0131 

dept of the a is force 


78 S00903DS 


AF WEAPONS LABOBATORi 

ATTN: TECHNICAL LIBRARY (SUL) 

KIRTLAND AFB NH 87117 


POSTMASTER; 






If Undeliverable (Section 158 
Postal Manual) l>o Not Return 






